ON HIGHER ORDER SUGAWARA OPERATORS 
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Abstract. The higher Sugawara operators acting on the Verma modules over the 
affine Kac-Moody algebra at the critical level are related to the higher Hamiltonians 
of the Gaudin model due to work of Feigin, Frenkel and Rcshctikhin. An explicit 
construction of the higher Hamiltonians in the case of gl n was given recently by 
Talalaev. We propose a new approach to these results from the viewpoint of the 
vertex algebra theory by proving directly the formulas for the higher order Sugawara 
operators. The eigenvalues of the operators in the Wakimoto modules of critical 
level are also calculated. 
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1. Introduction 

Let g be a simple Lie algebra over C. The corresponding affine Kac-Moody algebra 
is defined as the central extension 

(l.i) Q = Q[t,t- 1 ]®CK, flfor 1 ] = fl®c[t,r 1 ], 

where C [t, t~ l ] is the algebra of Laurent polynomials in t. Consider the Verma module 
M(A) over g at the critical level, so that the central element K acts as multiplication 
by the negative of the dual Coxeter number — h y . The Sugawara operators form a 
commuting family of g-endomorphisms of M(A). Such families of operators were first 
constructed by Goodman and Wallach [12] for the A series and independently by 
Hayashi [13] for the A, B, C series. These constructions were used in both papers for 
a derivation of the character formula for the irreducible quotient L(X) of M(A), under 
certain conditions on A, thus proving the Kac-Kazhdan conjecture. The existence 
of families of Sugawara operators for any simple Lie algebra g was established by 
Feigin and Frenkel [6] by providing a description of the center of the local completion 
U_ftv(g) loc of the universal enveloping algebra of g at the critical level. This work is 
based on the vertex algebra theory and makes use of the Wakimoto modules over g; 
see also [9]. 

On the other hand, due to work of Feigin, Frenkel and Reshetikhin [8] the central 
elements of the local completion U_/jv(g)i oc are related to the higher Hamiltonians of 
the Gaudin model describing quantum spin chain. In the case of g = g[ n , Talalaev [20] 
produced a remarkably simple explicit construction of the higher Hamiltonians. His 
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construction employs a close connection between the XXX and Gaudin quantum in- 
tegrable models which was extended further by Mukhin, Tarasov and Varchenko [18J. 
In [18] a version of the Bethe ansatz was used to produce some families of common 
eigenvectors of commutative subalgebras of the Yangian Y(g[ n ) and the universal 
enveloping algebra U(g[ n [t]) in tensor products of irreducible gl n -modules (for an 
introduction to the Yangians see [T7]). 

The relationship between the higher Gaudin Hamiltonians and the local completion 
of the universal enveloping algebra was employed in [3] to produce central elements of 
U-fr v (0( n )ioc in an explicit form. However, the proof given there is indirect and relies 
on [6], [IH] and [20] . In this paper we derive this result directly with the use of the 
vertex algebra theory. We construct two complete sets of higher Sugawara operators 
for g[„ which thus provides explicit formulas for the singular vectors of the Verma 
modules M(A) at the critical level. The first set is directly related to [20J and the 
construction recovers the corresponding results of [2], [3], as well as the main result 
of [20]. We also reproduce a description of the classical W-algebra associated to gl n 
and calculate the eigenvalues of the central elements in the Wakimoto modules. 

2. Preliminaries on vertex algebras 

2.1. Universal affme vertex algebra. A vertex algebra V is a vector space with 
the additional data (Y,T, 1), where the state-field correspondence Y is a map 

Y : V -> EndV^s,* -1 ]], 

the infinitesimal translation T is an operator T : V — > V, and 1 is a vacuum vector 
1 G V. These data must satisfy certain axioms; see e.g., [I], [10], [15]. For a G V we 
write 

Y(a, z) = } j a^z'™' 1 , a (n) G End V. 

In particular, for all a, b G V we must have a( n ) • b = for n ^> 0. The span in End V 
of all Fourier coefficients a( n ) of all vertex operators Y(a,z) is a Lie subalgebra of 
End V which we denote by U. The commutator of Fourier coefficients is given by 

(2' 1 ) [0(m),M = ( Z) ( a(n) ' h \rn + k-nY 

Note that (12.11) is equivalent to the formula for the commutator [Y(a, z),Y(b,w)] 
implied by the operator product expansion formula; see e.g. [1], pUl Chap. 3]. The 
center of the vertex algebra V is its commutative vertex subalgebra spanned by all 
vectors b G V such that a(„) - 6 = for all a G V and n ^ 0. The following observation 
will play an important role below. If b is an element of the center of V, then (12.11) 
implies [Y(a, z),Y(b,w)} = for all a G V. In other words, all Fourier coefficients 
belong to the center of the Lie algebra IA. 
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In this paper we are interested in the vertex algebras associated with affine Kac- 
Moody algebras g defined in (II. ip . where g is a finite-dimensional simple (or reductive) 
Lie algebra over C. Following [4], define the universal affine vertex algebra V(g) as 
the quotient of the universal enveloping algebra U(g) by the left ideal generated by 
g[t\. The vacuum vector is 1 and the translation operator is determined by 

(2.2) T:1^0, [T,K]=0 and [T, X[n\] = -nX[n - 1], leg, 

where we write X[n] = Xt n . The state-field correspondence Y is defined by setting 
r(M) = id, 

Y(K,z) = K, Y(J a [-l],z) = J a (z) = J a [r\z- r -\ 

where J 1 , . . . , J d is a basis of g, and then extending the map to the whole of V(g) 
with the use of normal ordering. Namely, the normally ordered product of fields 

a(z) = ^^a( r )2; _r_1 and b(w) = b^w~ r ~ 1 

rGZ r€Z 

is the formal power series 

(2.3) : a(z)b(w): = a{z)+b{w) + b(w)a(z)-, 
where 

a(z)+ = ^^a (r )Z~ r ~ 1 and a ( z )- — a( r )Z~ r-1 . 

r<0 r>0 

This definition extends to an arbitrary number of fields with the convention that the 
normal ordering is read from right to left. Then 

Y(J^[ ri }...J a -[r m },z) 

= 7 tt; ~7 ttt : d-^-\F^z) . . . d^-\F-{z) : 

(-ri - 1)! . . . (-r m - 1)! 

with r\ ^ ■ • • ^ r m < 0, and cij ^ a^+i if = r i+ i. 

The Lie algebra U(g)i oc spanned by the Fourier coefficients of the fields Y(a, z) 
with a G V(g) is known as the local completion of the universal enveloping algebra 
U(g) see [101 Sec. 3.5]. More precisely, the Lie algebra U(g)i oc is embedded into 
the completion U(g) of U(g) in the topology whose basis of the open neighborhoods 
of is formed by the left ideals generated by g ® t N< C\t\ with N ^ 

2.2. Segal— Sugawara vectors. As a vector space, V(g) can be identified with the 
space of polynomials in K with coefficients in the universal enveloping algebra U(g_), 
where g_ = g®t _1 C[t -1 ]. An element S G U(g_) C V(g) is called a Segal- Sugawara 
vector if 

(2.4) g[t]Se(K + h v )V(g), 
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where h v is the dual Coxeter number for g. A well-known example of such a vector 
is provided by the quadratic element 

(2-5) S = ~^J a [-l]J a [-l], 

0=1 

where J\, . . . , Jd is the basis of g dual to J 1 , . . . , J d with respect to the suitably 
normalized invariant bilinear form; see e.g. PHI Chap. 2]. 

The Segal-Sugawara vectors form a vector subspace of U(g_) which we will denote 
by Moreover, 3(g) is clearly closed under the multiplication in U(g_). The 

structure of the algebra 3(g) was described in [5J; see also [9J. In order to formulate 
the result, for any element S G U(g_) denote by S its image in the associated graded 
algebra grU(g_) = S(g_). We call the set of elements 

Si,...,S n G U(g_), n = rkg 

a complete set of Segal-Sugawara vectors if each Si satisfies (12. 4p and the correspond- 
ing elements Si, . . . , S n coincide with the images of certain algebraically independent 
generators of the algebra of invariants S(g) s under the embedding S(g) S(g_) 
defined by the assignment J a h-» J a [— 1]. In accordance to [6J, 3(g) is the algebra of 
polynomials in infinitely many variables, 

(2.6) z(&)=C[T r S l \l = l,...,n, r > 0]. 

A different proof of this result was given in a more recent work 0, Theorem 9.6]. 

Equivalently, the algebra 3(g) can be viewed as the center of the vertex algebra 
VCftv(g), where for each k G C the universal affine vertex algebra V k (q) of level k is 
the quotient of V(g) by the ideal (K — k)V(q). The space V K (g) is also called the 
vacuum representation of g of level k. The center of V K (g) is trivial (i.e., coincides 
with the multiples of the vacuum vector) for all values of k except for the critical 
value k = — h y . If k 7^ — h v then the vertex algebra V k (q) has the standard conformal 
structure determined by the Segal-Sugawara vector (12. 5ft . Due to the celebrated 
Sugawara construction, the Fourier coefficients of the field 

v ' a=l 

generate an action of the Virasoro algebra on the space V k (q). On the other hand, if 
k = —h v , then the Fourier coefficients of the field 

d 
a=l 

are the Sugawara operators, they commute with the action of g on the Verma mod- 
ule M(A) at the critical level. The existence of higher order Sugawara operators 
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on the Verma modules was established in p2] and [13], and certain complete (or 
fundamental) sets of such operators were produced for types A, B, C; see also pQ. 

Given a complete set of Segal-Sugawara vectors Si G V(fl) with I = 1, . . . ,n, the 
corresponding fields 

(2.7) S l (z)=Y(S l ,z)=J2Si,(r ) z- r ~ 1 

will form a complete set of Sugawara fields. This terminology can be shown to be 
consistent with [12] and [13]. The (higher) Sugawara operators are the Fourier coeffi- 
cients Si t ( r ) G U(fl)ioc of the fields (12.71) . If the highest weight A of the Verma module 
M(A) over g is of critical level, then K + h v is the zero operator on M(A). Due to (12.11) 
and (12 .4p . the Sugawara operators form a commuting family of g-endomorphisms of 
M(A). Moreover, an explicit construction of a complete set of Sugawara fields (12.71) 
yields an explicit description of all singular vectors of the generic Verma modules 
M(A) as polynomials in the Sugawara operators SW r ) with r < applied to the 
highest vector of M(A); see [131 Proposition 4.10] for a precise statement. 

The local completion at the critical level is the quotient U -hy (fl)ioc of U(g)i oc by the 
ideal generated by K + h y . The Sugawara operators can be viewed as the elements 
of the center of U -^v (g)i oc . In accordance to the description of the center which was 
originally given in [6], all central elements are obtained from 3(g) by the state-field 
correspondence. 

3. Segal-Sugawara vectors for gl n 

We let eij denote the standard basis elements of the Lie algebra gl n which we equip 
with the invariant symmetric bilinear form 

(X, Y) =ti(XY) - -trXtrV, X,Y e gl n . 

n 

The element en + ■ ■ • + e nn spans the kernel of the form, therefore it defines a non- 
degenerate invariant symmetric bilinear form on the Lie algebra sl n . The elements 
Cijt r of gl n [i, with r G Z will be denoted by ey[r]. Then the affine Kac-Moody 
algebra gt n = gl n [t,t _1 ] © CK has the commutation relations 

5 5 

(3.1) [etf[r],e w [s]] = 5 kj eu[r + s] - 5 u e kj [r + s] + K(5 kj 5 u - JL ^) rd r _ s , 

and the element K is central. The normalization of the form is chosen in such a way 
that the critical level — n coincides with the negative of the dual Coxeter number for 
sl n . We will also need the extended Lie algebra g[„ ©Cr, where for the element r we 
have the relations 



(3.2) 



[t, [r] ] = -r e {j [r - 1] , [r, K] = 0. 
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We identify the universal affine vertex algebra V(gl n ) with the vector space of 
polynomials in K with coefficients in U(g_) with g_ = t~ $ n \t~ ]■ We will also 
consider the quotient of U(g[ n © Cr) by the left ideal generated by fll n [£]. As a vector 
space, this quotient will be identified with V(Ql n ) © C[r], which may be viewed as 
the space of polynomials in r with coefficients in V(gl n ). 

For an arbitrary n x n matrix A = [a^] with entries in a ring we define its column- 
determinant cdet A by the formula 

(3.3) cdet A = ^ sgn a ■ a a{1)1 . . . a a{n)n , 

17 

summed over all permutations a of the set {1, . . . , n}. 

Consider the n x n matrix r + E[— 1] with entries in U(g_ 



r + E[-l] 



V + eu[-l] 
e 2 i[-l] 



ei 2 [-l] 
r + e 22 [-l] 



ein[ 
e 2n [ 



Cr) given by 
-1] 



e n i[-l] e n2 [-l] 



r + e nn [-l] 



Its column-determinant cdet(r + E[— 1]) is an element of U(g_ 
regard as an element of V(Ql n ) ® C[t]. 

Our main result is a direct proof of the following theorem. 



Cr) which we also 



Sn 



Theorem 3.1. The coefficients S±, . . . ,S n of the polynomial 
(3.4) cdetfr + E[-l}) = r n + S^- 1 + ■■■ + S n ^r 

form a complete set of Segal-Sugawara vectors in V(gi n ). Hence, 3(fll n ) is the algebra 
of polynomials, 



3 (0lr, 



'[T r Si | 1 = 1,..., n; 0]. 



We will prove the theorem in Section HI Here we consider some examples and 
corollaries, and derive the second construction of a complete set of Segal-Sugawara 
vectors; see Theorem 13.31 below. 

First we point out that regarding the Lie algebra sl n as the quotient of gl n by the 
relation en + • — h e nn = 0, we obtain the respective complete set of Segal-Sugawara 
vectors in V(sl n ). In particular, the vector S% vanishes, while — 5* 2 coincides with the 
vector (1231) . 

In what follows, the usual vertical line notation for determinants will always be 
used for column-determinants. 

Example 3.2. If n = 2 then 
S'i = e 11 [-l]+e 22 [-l], 

e n [-l] e 12 [-l] 
e 2 i[-l] e 22 [-l] 



S 9 



+ e 22 [-2] = en[-l] e 22 [-l] - e 21 [-l] e 12 [-l] + e 22 [-2]. 
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If n 



3 then 

S , i = e 11 [-l] + e 22 [-l]+e 3 3[-l], 



eil [-l] e 12 [-l] 
e 2 i[-l] e 22 [-l] 
+ e 22 [-2]+2e 33 [-2] 



eil [-l] e 13 [-l] 
e 3 i[-l] e 33 [-l] 



+ 



e 22 [-l] e 23 [-l] 
e 3 2[-l] e 33 [-l] 



and 



eu[-l] ei 2 [-l] e 13 [-l] 
e 2 i[-l] e 22 [-l] e 23 [-l] 
e 3 i[-l] e 32 [-l] e 33 [-l] 



2e 33 [-3] 



+ 



eni-1] e 13 [-2] 
e 3 i[-l] e 33 [-2] 



+ 



e 22 [-l] e 23 [-2] 
e 32 [-l] e 33 [-2] 



e 22 [-2] e 23 [-l] 
e 32 [-2] e 33 [-l] 



□ 



Now we use Theorem 13.11 to obtain the second complete set of Segal-Sugawara 
vectors in V(gl n ). 



T nn form a 



Theorem 3.3. For any k^ all coefficients T k i in the expansion 

tr( T + E[-l]) k = T k0 r k + T kl r k ~ x + ■ • • + T kk 

are Segal-Sugawara vectors in V(Ql n ). Moreover, the elements Tu, . 
complete set of Segal-Sugawara vectors. Hence, 3(fl[ n ) is the algebra of polynomials, 

M n ) = C[T r T u | l = l,...,n; r > 0]. 

Proof. Observe that if we replace r by u + r, where u is a complex variable, then 
relations (13.21) will still hold. Therefore, the coefficients in the expansion of the 
polynomial cdet(w +r + E[— 1]) in the powers of u and r are Segal-Sugawara vectors 
in V(Ql n ). Hence, the first part of the theorem is implied by the identity 



(3.5) d u cdet(u + T + E[-1]) = cdet(w + r + E[-l]) ^(~ r 



^-^trfr- 



E[-1)Y 



k=0 



since it provides an expression for the vectors T k i in terms of the Si regarded as 
elements of U(jj_). The identity (I3.5P can be viewed as a noncommutative analogue of 
the Liouville formula (as well as the Newton identities). It follows from [2, Theorem 4], 
due to Lemma fl~3l below F1 

In order to prove the second part, note that the elements Tu, . . . , T nn coincide with 



the images of the traces of powers tr E of the matrix E = [e^] with k 



n 



under the embedding S(g[ n ) S(fl_) defined by t— > ey[— 1]. However, the elements 
trS, . . . ,ti E n are algebraically independent and generate the algebra of invariants 



^^We thank Pavel Pyatov for pointing out that Theorem 4 in [2] can also be proved with the use 
of i?-matrix arguments; cf. [TJ], [T7J Chap. 7]. 
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Example 3.4. 


We have 
















Tio 


= n, 


T n 


= trE[- 


-1] 










T20 


= n, 


T 2 i 


= 2trE 


-1], 


T22 — 


tiE[- 


I] 2 +tiE[-2] 




T30 


= n, 


T31 


= 3tiE 


-1], 


T32 — 


3trE[ 


-l] 2 + 3tr£[- 


-2] 


and 


















T33 


= trE[- 


1] 3 + 


-2tr£[- 


1]E[ 


-2] +trE 


[-2}E 


-1] +2tiE[- 


3]. 



□ 

Remark 3.5. Some other families of Segal-Sugawara vectors can be constructed by 
using properties of Manin matrices. In particular, the quantum MacMahon master 
theorem [TT] leads to a construction of permanent-type vectors; see also [2]. □ 

4. Proof of Theorem 13.11 

We start by establishing some general properties of column-determinants. First 
we note that the column-determinant cdet A of a matrix A over an arbitrary ring 
changes sign if two rows of the matrix A are swapped. In particular, cdet A = if A 
has two identical rows. 

The next lemma is immediate from the definition of the column-determinant. 

Lemma 4.1. Let A = [a^] be an arbitrary n x n matrix with entries in a ring and b 
an element of the ring. Then 

an ■ ■ ■ [b, an] ■ ■ ■ a ln 
a 2 i ■ • • [b, a 2 j] ■ • • a 2n 



[6, cdet A] = 

a n i ■ ■ ■ [b, a n i\ ■ ■ ■ a. 

Lemma 14.11 implies one more property of column-determinants. 

Lemma 4.2. Let A = [a^] be an arbitrary n x n matrix with entries in a ring and 
b an element of the ring. Replace column j of A by the column whose all entries are 
zero except for the i-th entry equal to b. Then the column- determinant of this matrix 
can be written as 



a n ■ ■ 


■ ■ 


0>ln 




a u 




• a ln 




Oil ' ' 


• b ■ 




On 




• a in b 


a nl " ' 


■ ■ 






O-nl 


















a xl ■ ■ 


■ [6, aifc] • 


a ln 








+ (-iy +j 


n 

E 


a 2 i • ■ 


■ [6, a 2 k\ ■ 


a 2n 








k 


=3+1 




■ [b, a nk ] ■ 














O-nl ' ' 


0"nn 
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where the first determinant on the right hand side is obtained by moving column j to 
the last position, while row i and column j in the determinants occurring in the sum 
are deleted and the commutators [b,aik] occur in column k — 1. □ 

Now we recall some properties of a class of matrices introduced by Manin [16]. 
Following [2] we call a matrix A = [a^] with entries in a ring a Manin matrix if 

(hj aki - a M ctij = a kj a a - an a kj for all possible i, j, k, I. 

Such matrices are also known in the literature as right-quantum matrices (with q — 1); 



see [TT]. It is straightforward to verify that the column-determinant of a (square) 
Manin matrix will change sign if two columns are swapped. In particular, if a Manin 
matrix has two identical columns, then its column-determinant is zero. The next 
observation will be useful in the calculations below. 

Lemma 4.3. The matrix t+E[— 1] with entries in U(g_©Cr) is a Manin matrix. □ 

We begin proving Theorem 13.11 by verifying that the elements S±, . . . , S n satisfy 
QZ21) . Since 

[ eii [0],r]=0 and [e nn [l], r k ] = k r k ~ l e nn [0], 
it will be sufficient to check that for all i,j 

(4.1) eii [0]cdet(r + E[-l]) = and e nn [l] cdet(r + E[-l]) e (K + n)V(gl n ) 

in the $j[ n -module V(gl n ) <8> C[t\. The first relation in ( 14. ip is analogous to the well- 
known property of the Capelli determinant cdet[8 i: j(u— z+l)+e^]; this is a polynomial 
in u with coefficients in the center of the universal enveloping algebra U(gLj. This 
property can be verified by a direct argument, and we argue in a similar way to prove 
the relation. By Lemma [4. 1[ the polynomial e^ fO] cdet(r + E[— 1]J equals 





-e y [-l] 




o 


(4.2) 


••• e«[-l] - e#[-l] ••• 


+ E 


••• e ik [-l] ••• 












-e nj [-l] 




o 



where the dots indicate the same entries as in the column-determinant of the matrix 
t + E[— 1], except for column i in the first determinant and column k in the fc-th term 
in the sum; the entries shown in the middle belong to row j of each determinant. 
The first determinant can be written as the difference of the determinants 





(4.3) 



r + e«[-l] 




eiihl] 



e»i[-l] 
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Now, if % = j then the second determinant here equals cdet(r + E[— 1]), while the 
first determinant together with the sum over k in (14.21) equals cdet(r + £J[— 1]) by the 
analogue of the j-th row expansion formula for column-determinants. Hence the first 
relation in ( 14. ip holds in this case. 

If i ^ j, then the second determinant in ( 14.31) is obtained from cdet(r + E[— 1]) 
by replacing column i with column j. Similarly, by the row expansion formula for 
column-determinants, the first determinant in fl4.3[) together with the sum over k 
in ( 14. 2p is obtained from cdet(r + E[— 1]) by replacing row j with row i. Both 
determinants are zero due to Lemma 14.31 and the properties of column-determinants. 

Now we turn to the second relation in (14. ip . By (13. ip and (13.21) we have 

[e nn [1] , r + en [- 1] ] = e nn [0] - K' and [e nn [1] , e ni [- 1] ] = e ni [0] , 
where i ^ n and we have put K' = K/n, while 

[e nn [l], e in [-l) ] = -e in [0) and [e nn [l], r + e nn [-l] ] = (n - 1) K' + e nn [0]. 
Therefore, by Lemma fl~Tl the polynomial e nn [l) cdet(r + E[— 1]) equals 



(4.4) 



n-1 

E 

i=i 





e nn [0]-K' 


e n i[0] 



-em[0] 

-e in [0) 

-e n _i jn [0] 
(n-l)^' + e nn [0] 



where the dots replace the entries of the determinant cdet(r + E[— 1]), except for 
column i in the i-th summand and the last column in the last summand. Since each 
element ej n [0] annihilates the vacuum vector of V(gl n ), the last summand can be 
written as (n — 1)K' |r + E[— 1]| , where by |r + E[— 1] .. we denote the column- 
determinant of the matrix obtained from r + E[— 1] by deleting row and column i. 
Using this notation, we can represent (14.41) in the form 



n-1 



(4-5) £ 



i=l 





e n n[0] 
e n »[0] 



n-1 



K > J2\r + E[-l]\.. + (n - l)K' \t + E[-l]\ 



i=l 
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Our next step is to prove that the following key relation holds for each value of 
%— 1, . . . , n — 1: 







(4.6) 



e„„[0] 
en«[0] 



\t + E[-1]\-\t + E[-1\\ 



The determinant on the left hand side equals the sum of two determinants obtained 
by replacing e ni [0] or e nn [0] by 0, respectively Now we apply Lemma 1431 to each of 
these determinants. For the first one we have 





o 







(4.7) 




= (-l) n_< 


e nn [0] 




o 








n-l 



k=i+l 







e n fc[-l] 



-eij-l] 

-e n -i,n[-l] 




where the determinants in the second line are of the size (n— 1) x (n — 1) with row and 
column i deleted. Now add and subtract the following (n — 1) x (n — 1) determinant 
with row and column i deleted 



-r - e r , 



-1] 



\t + E[-1]\.. ■ (-T 

L 1 I vi, nn \ 



where r + E\— ll .. denotes the column-determinant of the matrix obtained from 

I L J I n,nn 

t + E[— 1] by deleting rows and columns i and n. Since the element e nn [0] annihilates 
the vacuum vector of V(gi n ), the expression for the determinant in f)4.7p simplifies to 



n-l 

E 

k=i+l 



\T + E[-l]\ ii + \T + E[-l}\(-T-e nn [-l})- 



e n fc[-l] 
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Applying now Lemma [4.21 in a similar way to the determinant in ( 14. 6 h with e nn [0] 
replaced with 0, we get 





o 







(4.8) 


o 


= (-I)""* 







•■■ e ni [0] ■■■ 




e n i[0) 



n-l 

+ (-ir* £ 


o 




-eii[-l] 


•■■ e nfe [-l] ■•• 


+ (-l) n ~ l 


C rari ,[ 1] 1] 


k=i+l 










o 




-e n _ij[-l] 



where, as before, the dots indicate the same entries as in the column-determinant 
cdet(r + E[— the determinants in the second line are of the size (n — 1) x (n— 1) 
with row n and column i deleted, and the entry e n f.[— 1] in the summation term occur 
in row % and column k — 1. Since the element e n j[0] annihilates the vacuum vector 
of V(gl n ), the first determinant on the right hand side of ( 14. 8 ft vanishes. In order 
to transform the remaining combination, add and subtract the (n — 1) x (n — 1) 
column-determinant 





-r 





obtained from cdet(r + E[— 1]) by deleting row n and column i and replacing the last 
column as indicated, with — r in the row i. After combining this difference with the 
last term in the expansion ( 14.81) we obtain the difference of two (n — 1) x (n — 1) 
column-determinants multiplied by (— l) n- \ where one of them is obtained from 
cdet(r + E[— 1]) by deleting row n and moving column % to replace the last column. 
By Lemma H~3l the corresponding (n— 1) x (n — 1) matrix is a Manin matrix. Hence, 
taking the signs into account we conclude that the determinant in question equals 
the minor |r+-B[— 1] For the remaining determinants we use the general property 
of column-determinants allowing us to permute rows. Moving row i in each of the 
determinants down to the last position and taking signs into account, we find that 
the determinant in (14.81) equals 



n-l 

E 

k=i+l 







e n fe[-l] 



E[-l]\ 



\t + E[-1]\.. (-t 
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Combining this with the expression for the determinant ( 14.7ft derived above, we com- 
plete the verification of ( 14.61) . 

Recalling now the expression ( 14.51) . we arrive at the relation 



n-l 



e nn [l] cdet(r + E[-l}) = — ^ ((n - 1) |r + E[-l]\ m - £ |r + 

i=i 

thus completing the proof of (14.11) . 

Finally, under the embedding S($j[ n ) <^-> S(g_) defined by i— > ey[— 1], the ele- 
ments Si, . . . ,S n coincide with the images of the respective coefficients of the charac- 
teristic polynomial det^-u + which are algebraically independent and generate 
the algebra of invariants S(gl n ) s[n . Therefore, Si, . . . , S n form a complete set of Segal- 
Sugawara vectors in V(gl n ). 



5. Center of the local completion and commutative subalgebras 



Applying Theorems 13.11 and 13.31 we can get complete sets of Sugawara fields in an 
explicit form. Consider the vertex algebra V(gl n ) and set t%j{z) = Y(eij[—l],z) so 
that 

-r-l 



eij{z) 



e ijl r \ z 



i,j = l,...,n. 



Introduce the n x n matrix d z + E(z) by 



d z + E(z) 



d z + en(z) ei 2 {z) ... e ln (z) 
e2i(z) d z + e 22 (z) ... e 2n (z) 



e n i(z) e n2 (z) ... d z + e nn (z)_ 

and expand its normally ordered column-determinant 

(5.1) : cdet(<9, + E{z)) : = d n z + S x {z) d n z ~ x + ■■■ + S n _i(z) d z + S n (z). 

Equivalently, the fields Si(z) are given by Si(z) = Y(Si, z), where the elements Si are 
defined in (ET3j) . 



Example 5.1. For n = 2 we have 
: cdet 



d z + en(z) e X2 (z) 
e 2 i(z) d z + e 22 (z) 



= ■ (d z + eu(z)) (d z + e 22 (z)) - e 2 i(z) e X2 {z) : 
so that 

Si(z) = eu(z) + e 22 (z), S 2 (z) = : e n (z) e 22 (z) - e 2X (z) ei 2 (z) : + e' 22 (z). 



□ 
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Introduce the fields T k i(z) = Y(T k i, z) corresponding to the Segal-Sugawara vectors 
Tki- More explicitly, they can be defined by the expansion of the normally ordered 
trace 

(5.2) : tr(d z + E(z)) k : = T k0 (z) d k z + T kl (z) d k ~ l + ■■■ + T kk (z) } k^O. 

Note that the Sugawara operators associated with the vectors T k i exhibit some 
similarity with the families of operators constructed in [12 J and [13], although the 
exact relation with those families is unclear. 

Recall that 3(flt n ) is the algebra of Segal-Sugawara vectors in V(gl n ) which can 
also be regarded as the center of the vertex algebra V- n (Ql n ). The center of the local 
completion U_ n (g[ n )i oc is the vector subspace 3(flt n ) which consists of the elements 
commuting with the action of gl n . It was proved in |6J that any element of 3(fO 
is a Fourier component of a field corresponding to an element of 3(gt„). By Theo- 
rems [3TT] and [33], the fields Y(S, z) with S G 3(gt„) can be interpreted as differential 
polynomials either in the fields Si(z), . . . , S n (z) or in the fields T u (z), . . . , T nn {z) 
with normally ordered products; see (15.11) and (15.21) . Hence, we obtain two explicit 
descriptions of 3(g( n ) formulated below. The first of them was originally given in [3], 
and the arguments there rely on the results of [6J, [IH] and [2D] . 

Corollary 5.2. The center 3(g(„) of the local completion U_„(gl n )i oc consists of 
the Fourier coefficients of all differential polynomials in either family of the fields 
S x {z),...,S n {z) orT u (z),. . . ,T nn (z). □ 

By the vacuum axiom of a vertex algebra, the application of the fields Si(z) and 
T k i(z) to the vacuum vector 1 of V- n ($l n ) yields power series in z which we denote 
respectively by 

(5.3) S l (z) + = J2 S «r) Z ~ r ~ 1 and T *l(z) + = J2 T t(rf~ r ~ 1 - 

r<0 r<0 

Since the Fourier coefficients of the fields Si(z) and T k i(z) belong to 3(flt n ), all coef- 
ficients of the series (15.31) belong to the center 3(fll n ) of the vertex algebra V- n (gl n ). 
By (15. ip and (15. 2p . the series (15.31) can be written in an explicit form with the use of 
the matrix E(z) + = [ey(z) + ]. We have 

cdet(d z + E(z)+) = d: + ^(z) + dr l + ■■■ + S n ^(z) + d z + S n {z) + 

and 

ti(d z + E(z) + ) k = T k0 (z) + d k z + T kl (z) + dt 1 + ■■■ + T kk (z)+. 
We arrive at the following result, whose first part dealing with the commuting family 
of the elements goes back to the original work [2D] ; see also [2] , [3] . 

Corollary 5.3. The elements of each of the families 
S^,, with I = 1, . . . ,n and r < 0, Tl, with ^ I ^ k and r < 0, 



SUGAWARA OPERATORS 



15 



belong to 3(f|l n ). Moreover, 3(gt„) is the algebra of polynomials 

3(2n) = C[5+ r) |/ = l,...,n, r<0]=C[T+ (r) | / = 1, . . . , n; r < 0]. 

Proof. We only need to prove the algebraic independence of the families of the ele- 
ments iSj ( r ) and TjJ,y This follows by comparing their highest degree components 
in the graded algebra grU(jj_) = S(g_) with those of the elements T r Si and T r Tu, 
respectively; see Theorems 13.11 and 13.31 In the general case such relationship between 
two families of generators of S(g_) was pointed out in [19]. □ 

Note that in our approach, the pairwise commutativity of the elements Sfr^ in 
U(g_) is a consequence of the fact that the are Segal-Sugawara vectors. The 
higher Gaudin Hamiltonians can be obtained by taking the images of the elements of 
the commutative subalgebra 3(gf n ) of U(g_) in the algebras U(Ql n )® m under certain 
evaluation homomorphisms; see e.g. [8], [20] for details. 

6. Eigenvalues in the Wakimoto modules 

As another application of the explicit formulas of Theorems 13.11 and 13.31 we recover 
a description of the classical W-algebra W(fll n ). This allows one to calculate the 
eigenvalues of the elements of the center 3(flt n ) m the Wakimoto modules of critical 
level. We start by recalling the main steps in the construction of these modules 
following [9]. 

Consider the Weyl algebra A(Ql n ) generated by the elements aij[r] with rGZ, 
i, j = 1, . . . , n and i ^ j and the defining relations 

l<kj[r], o«[s]] = hj 8 uK-s for * < 3, 

whereas all other pairs of the generators commute. The Fock representation M($j[„) 
of A(Qi n ) is generated by a vector |0) such that for i < j we have 

aij [r}\0) = 0, r^O and a^r] |0) = 0, r > 0. 

The vector space M(gl n ) carries a vertex algebra structure. In particular, |0) is the 
vacuum vector, and for i < j we have 

Y{ aij [-1] \0),z) = Y,^]^" 1 and Y(a 3i [0} \0),z) = aji [r] z~ r . 

We will denote these series by a^{z) and a^{z\ respectively. 
Take an n-tuple 

x(t) = {xi(t), X n(t)), xM = YxiMt-^ 1 e C((t)), 

where C((t)) denotes the algebra of formal Laurent series in t containing only a finite 
number of negative powers of t. The following formulas define a jj[ n -module structure 
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on the vector space M(gl n ): 

e iti+ i(z) h-> a iji+1 (z) + 2j : Pki a ki{z) ■ 

k<l 

(6.1) e ii( 2 ) ^ ^2 d L '■ aik(z)a H (z) : 

k<l 

e%+i,i(z) ^ 5^ : Qli a ki{z) : +Ci9 2 a i+lii (z) + - Xi+i{z)) a i+1>i (z), 

k<l 

where P k \ and Q l kl are certain polynomial expressions in the a pq (z) with p > q, while 
and c$ are certain coefficients. The precise formulas are obtained by applying the 
formulas [9, (4.6)-(4.8)] to the Lie algebra sl n and extending the action to gl n by 
setting 

eu(z) H h e nn (z) i-> xi{ z ) H ^ Xn(z)- 

These are the Wakimoto modules of critical level, denoted by W x m. 

Example 6.1. For n = 2 the explicit formulas are 

e 12 (z) h-> 012(2;) 

en(z) h-> - : 021(2;) 012(2;) : + 
e 22 (2;) 1 : 021(2;) 012(2;) : + ^2(2;) 

^21 (2;) ' - : a 2 i(2;) 2 012(2;) : - 2 d z a 2 i(z) + (xi(z) - X2(z)) a 2 i(z). 

The key fact leading to the construction of the Wakimoto modules is the existence 
of the vertex algebra homomorphism 

(6.3) p:V_ n (Ql n )^M(Ql n )®7r Q , 

where ttq is the algebra of polynomials 

7r = C[bi[r] \i = l,...,n; r < 0] 

in the variables Oj[r] which we consider as a commutative vertex algebra. The trans- 
lation operator on 7r is defined by 

T1 = 0, [T,bi[r)]=-rbi[r-l). 

The homomorphism (I6.3P is defined by the formulas (16. ip . where the series Xi( z ) 
should be respectively replaced by 

b l (z) = J2b l [r}z-r-\ 

r<0 

The image of the center $(&l n ) of the vertex algebra V_ n (gl n ) under the homomor- 
phism p is contained in ttq = 1 <8>7ro. Moreover, this image coincides with the classical 
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W -algebra W(gLj; see P, Sec. 9]. In order to recall the definition of W(gl n ), intro- 
duce the operators 

Qi : 7r -> 7r , i = l, ...,n-l, 

by the formulas 

bi(A) / 9 9 



ft = EE 



, A _ V dh [-r - 1] <96 l+1 [-r - 1] 

where the following notation was used. The second sum is taken over all partitions 
A = (Ai, . . . , A p ) of r, so that Ai ^ • • • ^ \ p > and Xx + ■ • • + X p — r. Furthermore, 

MA) = (6i[-Ai] - 6i4-i[-Ai]) . . . (bi[-Xp] - b i+1 [-X p }) 

and 

z x = l mi m 1 !2 m2 m 2 !...r m 'm r !, 

where denotes the multiplicity of k in A for each k. By definition, W(gl n ) consists 
of the elements of ttq, annihilated by all operators Qi, 

W(gl n )= p| KerQ,. 

1^ i ^n— 1 

Applying Theorems 13.11 and 13.31 we recover an explicit description of the algebra 
W(g[„); see e.g. [5], [Tj. In order to formulate the result, consider the extended 
algebra, which is isomorphic to tcq ® C [ vector space, with the relations 

[r,6i[r]] = -rbi[r - I). 

Corollary 6.2. The image of the column- determinant cdet(r + E[— 1]) under the 
homomorphism f !6.3p is given by 

p : cdet(r + E[-l]) h-> (r + 6 n [-l]) ■ • ■ (r + 6i[-l]). 

Hence, if the elements B-i are defined by 

(r + 6 n [-l]) • ■ ■ (r + b x [-l\) =r n + B 1 r"- 1 + • • • + B n , 

then W(gt n ) is the algebra of polynomials in the variables T r B-i with i = 1, . . . , n and 
r ^ 0. 

Proof. Using Lemma 14.31 and the properties of the Manin matrices, write the rows 
and columns of the matrix r + E[— 1] in their reverse orders so that the column 
determinant will take the form 

(6.4) cdet(r + £[-l]) = J^sgna ■ (S a{n)n r + e a(n)n [-l}) . . . (5 a{1)1 r + e a{1)1 [-l]) . 

a 

Now we apply the explicit description of the homomorphism (16.31) given in [9] . Given 
that the image of 3(gl n ) under this homomorphism is contained in 7To, we only need 
to keep the terms containing the elements bi[r\. The images of the generators ey[r] 
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with i < j and r < do not depend on the bi[r] and the images of the summands 
in ( 16. 4j) with <j(n) < n do not contribute to the image of the column determinant. 
Hence, the nonzero contribution is only provided by the summands in ( 16. 4p with 
cr(n) = n. Considering now a{n — 1) and arguing by induction, we conclude that the 
only summand in f ]6.4[) providing a nonzero contribution to the image corresponds to 
the identity permutation a. This proves the first part. The second statement follows 
from [9, Theorem 9.6]. □ 

Note that an alternative calculation of the image of cdet(r + E[— 1]) under the ho- 
momorphism (16. 3p can be performed with the use of the original expression provided 
by ( 13.31) and the application of [9, Theorem 11.3]. 

Example 6.3. If n = 2 then 

B 1 = b 1 [-l]+b 2 [-l], 

B 2 = b 1 [-l]b 2 [-l]+b 1 [-2]. 

If n = 3 then 

B 1 = 6 1 [-l]+6 2 [-l]+6 3 [-l], 

B 2 = b 2 [-l] + b x [-l\ 6 3 [-l] + M-l] 6 3 [-l] + 2 ^[-2] + b 2 [-2), 

B 3 = hi-1] b 2 [-l] b 3 [-l] + b 1 [-2] b 2 [-l] 

+ M-2]M- 1] + M-i]M-2] + 2 M -3]. 

Applying the identity (13.51) and Corollary 16.21 we can also calculate the images of 
the Segal-Sugawara vectors T w . 

Corollary 6.4. The images of the elements tv(r+E[— l]) k under the homomorphism 
(16.31) are found from the formula 

p : tk tr (r + E[-l]) k ^ ^ (l - t (r + ^[-1])) " • ■ • (l - t (r + 6*[-l]))~ 

k=0 i=l 

x(l-t(r + fei-if-l])) • • • (l - t (r + 6i[-l])), 

where t is a complex variable. 

Proof. The formula follows from (13.51) by calculating the derivative on the left hand 
side and then replacing u by —t~ l . □ 

The elements of the center 3(fl[ n ) of U_ n (g[ n )i oc act on the Wakimoto modules 
W x (t) as multiplications by scalars which can be calculated by using Corollaries 16.21 
and 16.41 Recall the description of 3(fl[„) provided by Corollary 15.21 
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Corollary 6.5. The coefficients of : cdet(d z + E(z)) : and : tr(<9 2 + E(z)^ k : act on 
W x m as multiplications by scalars found from the respective formulas 

: cdet{8 z + E{z)) : i-> + X n(») • • • + 

and 



^t fc :tr(c} 2 + E(z)) fc : ^ £ (l - t (d z + X i (*))) ■ ■ • (l - * (9. + 

f(0* + Xi_i(z)))---(l-*(d* + Xi(z))). 



fe=0 i=l 

• (1 



Acknowledgments 

We acknowledge the support of the Australian Research Council. The first author 
is grateful to the University of Sydney for hospitality during his visit. His research 
is also supported by the grant of Support for the Scientific Schools 8004.2006.2, 
RFBR grant 08-02-00287a, the ANR grant GIMP (Geometry and Integrability in 
Mathematics and Physics) and the grant MK-5056. 2007.1. 

References 

[1] F. A. Bais, P. Bouwknegt, M. Surridge and K. Schoutens, Extensions of the Virasoro algebra 
constructed from Kac-Moody algebras using higher order Casimir invariants, Nucl. Phys. B304 
(1988), 348-370. 

[2] A. Chervov and G. Falqui, Manin matrices and Talalaev's formula, J. Phys. A: Math. Theor. 

41 (2008), 194006 (28pp). 
[3] A. Chervov and D. Talalaev, Quantum spectral curves, quantum integrable systems and the 

geometric Langlands correspondence, preprint arXiv:hep-th/0604128. 
[4] A. De Sole and V. G. Kac, Finite vs affine W '-algebras, Japan. J. Math., 1 (2006), 137-261. 
[5] L. A. Dickey, Lectures on classical W -algebras, Acta Appl. Math. 47 (1997), 243-321. 
[6] B. Feigin and E. Frenkel, Affine Kac-Moody algebras at the critical level and Gelfand-Dikii 

algebras, Int. J. Mod. Phys. A7, Suppl. 1A (1992), 197-215. 
[7] B. Feigin and E. Frenkel, Integrals of motion and quantum groups, in "Integrable systems and 

quantum groups" (Montecatini Terme, 1993), 349-418, Lecture Notes in Math., 1620, Springer, 

Berlin, 1996. 

[8] B. Feigin, E. Frenkel and N. Reshetikhin, Gaudin model, Bethe ansatz and critical level, Comm. 

Math. Phys. 166 (1994), 27-62. 
[9] E. Frenkel, Wakimoto modules, opers and the center at the critical level, Adv. Math. 195 (2005), 

297-404. 

[10] E. Frenkel and D. Ben-Zvi, Vertex algebras and algebraic curves, Mathematical Surveys and 
Monographs, 88. American Mathematical Society, Providence, RI, 2001. 

[11] S. Garoufalidis, Thang T. Q. Le and D. Zeilberger, The quantum MacMahon master theorem, 
Proc. Nat. Acad. Sci. USA 103 (2006), 13928-13931. 

[12] R. Goodman and N. Wallach, Higher-order Sugawara operators for affine Lie algebras, Trans. 
Amer. Math. Soc. 315 (1989), 1-55. 

[13] T. Hayashi, Sugawara operators and Kac-Kazhdan conjecture, Invent. Math. 94 (1988), 13-52. 



20 



A. V. CHERVOV AND A. I. MOLEV 



[14] A. Isaev, O. Ogievetsky and P. Pyatov, On quantum matrix algebras satisfying the Cayley- 

Hamilton- Newton identities, J. Phys. A 32 (1999), L115-L121. 
[15] V. Kac, Vertex algebras for beginners, University Lecture Series, 10. American Mathematical 

Society, Providence, RI, 1997. 
[16] Yu. I. Manin, Quantum groups and noncommutative geometry, Universite de Montreal, Centre 

de Recherches Mathematiques, Montreal, QC, 1988. 
[17] A. Molev, Yangians and classical Lie algebras, Mathematical Surveys and Monographs, 143. 

American Mathematical Society, Providence, RI, 2007. 
[18] E. Mukhin, V. Tarasov and A. Varchenko, Bethe eigenvectors of higher transfer matrices, J. 

Stat. Mech. Theory Exp. 2006, no. 8, P08002, 44 pp. 
[19] L. G. Rybnikov, Uniqueness of higher Gaudin hamiltonians, preprint arXiv:math/0608588. 
[20] D. V. Talalaev, The quantum Gaudin system, Funct. Anal. Appl. 40 (2006), 73-77. 

Institute for Theoretical and Experimental Physics, Moscow 117259, Russia 
E-mail address: chervov@itep.ru 

School of Mathematics and Statistics, University of Sydney, NSW 2006, Australia 
E-mail address: alexm@matlis.usyd.edu.au 



